Single-cell RNA sequencing of mid-to-late stage spider embryos: new insights into spider development.

BACKGROUND: The common house spider Parasteatoda tepidariorum represents an emerging new model organism of arthropod evolutionary and developmental (EvoDevo) studies. Recent technical advances have resulted in the first single-cell sequencing (SCS) data on this species allowing deeper insights to be gained into its early development, but mid-to-late stage embryos were not included in these pioneering studies. RESULTS: Therefore, we performed SCS on mid-to-late stage embryos of Parasteatoda and characterized resulting cell clusters by means of in-silico analysis (comparison of key markers of each cluster with previously published information on these genes). In-silico prediction of the nature of each cluster was then tested/verified by means of additional in-situ hybridization experiments with additional markers of each cluster. CONCLUSIONS: Our data show that SCS data reliably group cells with similar genetic fingerprints into more or less distinct clusters, and thus allows identification of developing cell types on a broader level, such as the distinction of ectodermal, mesodermal and endodermal cell lineages, as well as the identification of distinct developing tissues such as subtypes of nervous tissue cells, the developing heart, or the ventral sulcus (VS). In comparison with recent other SCS studies on the same species, our data represent later developmental stages, and thus provide insights into different stages of developing cell types and tissues such as differentiating neurons and the VS that are only present at these later stages.

New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore.

BACKGROUND: Early during onychophoran development and prior to the formation of the germ band, a posterior tissue thickening forms the posterior pit. Anterior to this thickening forms a groove, the embryonic slit, that marks the anterior-posterior orientation of the developing embryo. This slit is by some authors considered the blastopore, and thus the origin of the endoderm, while others argue that the posterior pit represents the blastopore. This controversy is of evolutionary significance because if the slit represents the blastopore, then this would support the amphistomy hypothesis that suggests that a slit-like blastopore in the bilaterian ancestor evolved into protostomy and deuterostomy. RESULTS: In this paper, we summarize our current knowledge about endoderm and mesoderm development in onychophorans and provide additional data on early endoderm- and mesoderm-determining marker genes such as Blimp, Mox, and the T-box genes. CONCLUSION: We come to the conclusion that the endoderm of onychophorans forms prior to the development of the embryonic slit, and thus that the slit is not the primary origin of the endoderm. It is thus unlikely that the embryonic slit represents the blastopore. We suggest instead that the posterior pit indeed represents the lips of the blastopore, and that the embryonic slit (and surrounding tissue) represents a morphologically superficial archenteron-like structure. We conclude further that both endoderm and mesoderm development are under control of conserved gene regulatory networks, and that many of the features found in arthropods including the model Drosophila melanogaster are likely derived.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication.

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Developmental and genomic insight into the origin of the tardigrade body plan.

Tardigrada is an ancient lineage of miniaturized animals. As an outgroup of the well-studied Arthropoda and Onychophora, studies of tardigrades hold the potential to reveal important insights into body plan evolution in Panarthropoda. Previous studies have revealed interesting facets of tardigrade development and genomics that suggest that a highly compact body plan is a derived condition of this lineage, rather than it representing an ancestral state of Panarthropoda. This conclusion was based on studies of several species from Eutardigrada. We review these studies and expand on them by analyzing the publicly available genome and transcriptome assemblies of Echiniscus testudo, a representative of Heterotardigrada. These new analyses allow us to phylogenetically reconstruct important features of genome evolution in Tardigrada. We use available data from tardigrades to interrogate several recent models of body plan evolution in Panarthropoda. Although anterior segments of panarthropods are highly diverse in terms of anatomy and development, both within individuals and between species, we conclude that a simple one-to-one alignment of anterior segments across Panarthropoda is the best available model of segmental homology. In addition to providing important insight into body plan diversification within Panarthropoda, we speculate that studies of tardigrades may reveal generalizable pathways to miniaturization.

Expression of posterior Hox genes and opisthosomal appendage development in a mygalomorph spider.

Spiders represent an evolutionary successful group of chelicerate arthropods. The body of spiders is subdivided into two regions (tagmata). The anterior tagma, the prosoma, bears the head appendages and four pairs of walking legs. The segments of the posterior tagma, the opisthosoma, either lost their appendages during the course of evolution or their appendages were substantially modified to fulfill new tasks such as reproduction, gas exchange, and silk production. Previous work has shown that the homeotic Hox genes are involved in shaping the posterior appendages of spiders. In this paper, we investigate the expression of the posterior Hox genes in a tarantula that possesses some key differences of posterior appendages compared to true spiders, such as the lack of the anterior pair of spinnerets and a second set of book lungs instead of trachea. Based on the observed differences in posterior Hox gene expression in true spiders and tarantulas, we argue that subtle changes in the Hox gene expression of the Hox genes abdA and AbdB are possibly responsible for at least some of the morphological differences seen in true spiders versus tarantulas.

Comment on "The lower Cambrian lobopodian Cardiodictyon resolves the origin of euarthropod brains".

Strausfeld et al. (Report, 24 Nov 2022, p. 905) claim that Cambrian fossilized nervous tissue supports the interpretation that the ancestral panarthropod brain was tripartite and unsegmented. We argue that this conclusion is unsupported, and developmental data from living onychophorans contradict it.

Expression of netrin and its receptors uncoordinated-5 and frazzled in arthropods and onychophorans suggests conserved and diverged functions in neuronal pathfinding and synaptogenesis.

BACKGROUND: Development of the nervous system and the correct connection of nerve cells require coordinated axonal pathfinding through an extracellular matrix. Outgrowing axons exhibit directional growth toward or away from external guidance cues such as Netrin. Guidance cues can be detected by growth cones that are located at the end of growing axons through membrane-bound receptors such as Uncoordianted-5 and Frazzled. Binding of Netrin causes reformation of the cytoskeleton and growth of the axon toward (or away from) the source of Netrin production. RESULTS: Here, we investigate the embryonic mRNA expression patterns of netrin genes and their potential receptors, uncoordinated-5 and frazzled in arthropod species that cover all main branches of Arthropoda, that is, Pancrustacea, Myriapoda, and Chelicerata. We also studied the expression patterns in a closely related outgroup species, the onychophoran Euperipatoides kanangrensis, and provide data on expression profiles of these genes in larval tissues of the fly Drosophila melanogaster including the brain and the imaginal disks. CONCLUSION: Our data reveal conserved and diverged aspects of neuronal guidance in Drosophila with respect to the other investigated species and suggest a conserved function in nervous system patterning of the developing appendages.

Review of extra-embryonic tissues in the closest arthropod relatives, onychophorans and tardigrades.

The so-called extra-embryonic tissues are important for embryonic development in many animals, although they are not considered to be part of the germ band or the embryo proper. They can serve a variety of functions, such as nutrient uptake and waste removal, protection of the embryo against mechanical stress, immune response and morphogenesis. In insects, a subgroup of arthropods, extra-embryonic tissues have been studied extensively and there is increasing evidence that they might contribute more to embryonic development than previously thought. In this review, we provide an assessment of the occurrence and possible functions of extra-embryonic tissues in the closest arthropod relatives, onychophorans (velvet worms) and tardigrades (water bears). While there is no evidence for their existence in tardigrades, these tissues show a remarkable diversity across the onychophoran subgroups. A comparison of extra-embryonic tissues of onychophorans to those of arthropods suggests shared functions in embryonic nutrition and morphogenesis. Apparent contribution to the final form of the embryo in onychophorans and at least some arthropods supports the hypothesis that extra-embryonic tissues are involved in organogenesis. In order to account for this role, the commonly used definition of these tissues as 'extra-embryonic' should be reconsidered. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns.

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

Embryonic expression patterns of Wnt genes in the RTA-clade spider Cupiennius salei.

Spiders represent widely used model organisms for chelicerate and even arthropod development and evolution. Wnt genes are important and evolutionary conserved factors that control and regulate numerous developmental processes. Recent studies comprehensively investigated the complement and expression of spider Wnt genes revealing conserved as well as diverged aspects of their expression and thus (likely) function among different groups of spiders representing Mygalomorphae (tarantulas), and both main groups of Araneae (true spiders) (Haplogynae/Synspermiata and Entelegynae). The allegedly most modern/derived group of entelegyne spiders is represented by the RTA-clade of which no comprehensive data on Wnt expression were available prior to this study. Here, we investigated the embryonic expression of all Wnt genes of the RTA-clade spider Cupiennius salei. We found that most of the Wnt expression patterns are conserved between Cupiennius and other spiders, especially more basally branching species. Surprisingly, most differences in Wnt gene expression are seen in the common model spider Parasteatoda tepidariorum (a non-RTA clade entelegyne species). These results show that data and conclusions drawn from research on one member of a group of animals (or any other organism) cannot necessarily be extrapolated to the group as a whole, and instead highlight the need for comprehensive taxon sampling.

Phylogenetic analysis of forkhead transcription factors in the Panarthropoda.

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.

Lack of evidence for conserved parasegmental grooves in arthropods.

In the arthropod model species Drosophila melanogaster, a dipteran fly, segmentation of the anterior-posterior body axis is under control of a hierarchic gene cascade. Segmental boundaries that form morphological grooves are established posteriorly within the segmental expression domain of the segment-polarity gene (SPG) engrailed (en). More important for the development of the fly, however, are the parasegmental boundaries that are established at the interface of en expressing cells and anteriorly adjacent wingless (wg) expressing cells. In Drosophila, both segmental and transient parasegmental grooves form. The latter are positioned anterior to the expression of en. Although the function of the SPGs in establishing and maintaining segmental and parasegmental boundaries is highly conserved among arthropods, parasegmental grooves have only been reported for Drosophila, and a spider (Cupiennius salei). Here, we present new data on en expression, and re-evaluate published data, from four distantly related spiders, including Cupiennius, and a distantly related chelicerate, the harvestman Phalangium opilio. Gene expression analysis of en genes in these animals does not corroborate the presence of parasegmental grooves. Consequently, our data question the general presence of parasegmental grooves in arthropods.

Widespread retention of ohnologs in key developmental gene families following whole-genome duplication in arachnopulmonates.

Whole-genome duplications (WGDs) have occurred multiple times during animal evolution, including in lineages leading to vertebrates, teleosts, horseshoe crabs, and arachnopulmonates. These dramatic events initially produce a wealth of new genetic material, generally followed by extensive gene loss. It appears, however, that developmental genes such as homeobox genes, signaling pathway components and microRNAs are frequently retained as duplicates (so-called ohnologs) following WGD. These not only provide the best evidence for WGD, but an opportunity to study its evolutionary consequences. Although these genes are well studied in the context of vertebrate WGD, similar comparisons across the extant arachnopulmonate orders are patchy. We sequenced embryonic transcriptomes from two spider species and two amblypygid species and surveyed three important gene families, Hox, Wnt, and frizzled, across these and 12 existing transcriptomic and genomic resources for chelicerates. We report extensive retention of putative ohnologs, further supporting the ancestral arachnopulmonate WGD. We also found evidence of consistent evolutionary trajectories in Hox and Wnt gene repertoires across three of the six arachnopulmonate orders, with interorder variation in the retention of specific paralogs. We identified variation between major clades in spiders and are better able to reconstruct the chronology of gene duplications and losses in spiders, amblypygids, and scorpions. These insights shed light on the evolution of the developmental toolkit in arachnopulmonates, highlight the importance of the comparative approach within lineages, and provide substantial new transcriptomic data for future study.

A chelicerate Wnt gene expression atlas: novel insights into the complexity of arthropod Wnt-patterning.

The Wnt genes represent a large family of secreted glycoprotein ligands that date back to early animal evolution. Multiple duplication events generated a set of 13 Wnt families of which 12 are preserved in protostomes. Embryonic Wnt expression patterns (Wnt-patterning) are complex, representing the plentitude of functions these genes play during development. Here, we comprehensively investigated the embryonic expression patterns of Wnt genes from three species of spiders covering both main groups of true spiders, Haplogynae and Entelegynae, a mygalomorph species (tarantula), as well as a distantly related chelicerate outgroup species, the harvestman Phalangium opilio. All spiders possess the same ten classes of Wnt genes, but retained partially different sets of duplicated Wnt genes after whole genome duplication, some of which representing impressive examples of sub- and neo-functionalization. The harvestman, however, possesses a more complete set of 11 Wnt genes but with no duplicates. Our comprehensive data-analysis suggests a high degree of complexity and evolutionary flexibility of Wnt-patterning likely providing a firm network of mutational protection. We discuss the new data on Wnt gene expression in terms of their potential function in segmentation, posterior elongation, and appendage development and critically review previous research on these topics. We conclude that earlier research may have suffered from the absence of comprehensive gene expression data leading to partial misconceptions about the roles of Wnt genes in development and evolution.

Panarthropod tiptop/teashirt and spalt orthologs and their potential role as "trunk"-selector genes.

BACKGROUND: In the vinegar fly Drosophila melanogaster, the homeodomain containing transcription factor Teashirt (Tsh) appears to specify trunk identity in concert with the function of the Hox genes. While in Drosophila there is a second gene closely related to tsh, called tiptop (tio), in other arthropods species only one copy exists (called tio/tsh). The expression of tsh and tio/tsh, respectively, is surprisingly similar among arthropods suggesting that its function as trunk selector gene may be conserved. Other research, for example on the beetle Tribolium castaneum, questions even conservation of Tsh function among insects. The zinc-finger transcription factor Spalt (Sal) is involved in the regulation of Drosophila tsh, but this regulatory interaction does not appear to be conserved in Tribolium either. Whether the function and interaction of tsh and sal as potential trunk-specifiers, however, is conserved is still unclear because comparative studies on sal expression (except for Tribolium) are lacking, and functional data are (if at all existing) restricted to Insecta. RESULTS: Here, we provide additional data on arthropod tsh expression, show the first data on onychophoran tio/tsh expression, and provide a comprehensive investigation on sal expression patterns in arthropods and an onychophoran. CONCLUSIONS: Our data support the idea that tio/tsh genes are involved in the development of "trunk" segments by regulating limb development. Our data suggest further that the function of Sal is indeed unlikely to be conserved in trunk vs head development like in Drosophila, but early expression of sal is in line with a potential homeotic function, at least in Arthropoda.

Molecular evidence for a single origin of ultrafiltration-based excretory organs.

Excretion is an essential physiological process, carried out by all living organisms, regardless of their size or complexity.1-3 Both protostomes (e.g., flies and flatworms) and deuterostomes (e.g., humans and sea urchins) possess specialized excretory organs serving that purpose. Those organs exhibit an astonishing diversity, ranging from units composed of just few distinct cells (e.g., protonephridia) to complex structures, built by millions of cells of multiple types with divergent morphology and function (e.g., vertebrate kidneys).4,5 Although some molecular similarities between the development of kidneys of vertebrates and the regeneration of the protonephridia of flatworms have been reported,6,7 the molecular underpinnings of the development of excretory organs have never been systematically studied in a comparative context.4 Here, we show that a set of transcription factors (eya, six1/2, pou3, sall, lhx1/5, and osr) and structural proteins (nephrin, kirre, and zo1) is expressed in the excretory organs of a phoronid, brachiopod, annelid, onychophoran, priapulid, and hemichordate that represent major protostome lineages and non-vertebrate deuterostomes. We demonstrate that the molecular similarity observed in the vertebrate kidney and flatworm protonephridia6,7 is also seen in the developing excretory organs of those animals. Our results show that all types of ultrafiltration-based excretory organs are patterned by a conserved set of developmental genes, an observation that supports their homology. We propose that the last common ancestor of protostomes and deuterostomes already possessed an ultrafiltration-based organ that later gave rise to the vast diversity of extant excretory organs, including both proto- and metanephridia.

Oscillating waves of Fox, Cyclin and CDK gene expression indicate unique spatiotemporal control of cell cycling during nervous system development in onychophorans.

Forkhead box (Fox) genes code for a class of transcription factors with many different fundamental functions in animal development including cell cycle control. Other important factors of cell cycle control are Cyclins and Cyclin-dependent kinases (CDKs). Here we report on the oscillating expression of three Fox genes, FoxM, FoxN14 (jumeaux) and FoxN23 (Checkpoint suppressor like-1), Cyclins and CDKs in an onychophoran, a representative of a relatively small group of animals that are closely related to the arthropods. Expression of these genes is in the form of several waves that start as dot-like domains in the center of each segment and then transform into concentric rings that run towards the periphery of the segments. This oscillating gene expression, however, occurs exclusively along the anterior-posterior body axis in the tissue ventral to the base of the appendages, a region where the central nervous system and the enigmatic ventral and preventral organs of the onychophoran develop. We suggest that the oscillating gene expression and the resulting waves of expression we report are likely correlated with cell cycle control during the development of the onychophoran nervous system. This intriguing patterning appears to be unique for onychophorans as it is not found in any of the arthropods we also investigated in this study, and is likely correlated with the slow embryonic development of onychophorans compared to arthropods.

Expression of the zinc finger transcription factor Sp6-9 in the velvet worm Euperipatoides kanangrensis suggests a conserved role in appendage development in Panarthropoda.

The Sp-family genes encode important transcription factors in animal development. Here we investigate the embryonic expression patterns of the complete set of Sp-genes in the velvet worm Euperipatoides kanangrensis (Onychophora), with a special focus on the Sp6-9 ortholog. In arthropods, Sp6-9, the ortholog of the Drosophila melanogaster D-Sp1 gene plays a conserved role in appendage development. Our data show that the expression of Sp6-9 during the development of the velvet worm is conserved, suggesting that the key function of the Sp6-9 gene dates back to at least the last common ancestor of arthropods and onychophorans and thus likely the last common ancestor of Panarthropoda.

The embryonic expression pattern of a second, hitherto unrecognized, paralog of the pair-rule gene sloppy-paired in the beetle Tribolium castaneum.

In the fly Drosophila melanogaster, a hierarchic segmentation gene cascade patterns the anterior-posterior body axis of the developing embryo. Within this cascade, the pair-rule genes (PRGs) transform the more uniform patterning of the higher-level genes into a metameric pattern that first represents double-segmental units, and then, in a second step, represents a true segmental pattern. Within the PRG network, primary PRGs regulate secondary PRGs that are directly involved in the regulation of the next lower level, the segment-polarity genes (SPGs). While the complement of primary PRGs is different in Drosophila and the beetle Tribolium, another arthropod model organism, both paired (prd) and sloppy-paired (slp), acts as secondary PRGs. In earlier studies, the interaction of PRGs and the role of the single slp ortholog in Tribolium have been investigated in some detail revealing conserved and diverged aspects of PRG function. In this study, I present the identification and the analysis of embryonic expression patterns of a second slp gene (called slp2) in Tribolium. While the previously identified gene, slp, is expressed in a typical PRG pattern, expression of slp2 is more similar to that of the downstream-acting SPGs, and shows expression similarities to slp2 in Drosophila. The previously reported differences between the function of slp in Drosophila and Tribolium may partially account for the function of the newly identified second slp paralog in Tribolium, and it may therefore be advised to conduct further studies on PRG function in the beetle.

The forkhead box containing transcription factor FoxB is a potential component of dorsal-ventral body axis formation in the spider Parasteatoda tepidariorum.

In the spider, determination of the dorsal-ventral body (DV) axis depends on the interplay of the dorsal morphogen encoding gene decapentaplegic (Dpp) and its antagonist, short gastrulation (sog), a gene that is involved in the correct establishment of ventral tissues. Recent work demonstrated that the forkhead domain encoding gene FoxB is involved in dorsal-ventral axis formation in spider limbs. Here, Dpp likely acts as a dorsal morphogen, and FoxB is likely in control of ventral tissues as RNAi-mediated knockdown of FoxB causes dorsalization of the limbs. In this study, we present phenotypes of FoxB knockdown that demonstrate a function in the establishment of the DV body axis. Knockdown of FoxB function leads to embryos with partially duplicated median germ bands (Duplicitas media) that are possibly the result of ectopic activation of Dpp signalling. Another class of phenotypes is characterized by unnaturally slim (dorsal-ventrally compressed) germ bands in which ventral tissue is either not formed, or is specified incorrectly, likely a result of Dpp over-activity. These results suggest that FoxB functions as an antagonist of Dpp signalling during body axis patterning, similarly as it is the case in limb development. FoxB thus represents a general player in the establishment of dorsal-ventral structures during spider ontogeny.